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We investigate the time evolution of a superposition of 
macroscopically distinct quantum states in a system of two- 
level atoms interacting with a thermal environment of pho- 
ton modes. We show that the atomic coherent states are 
robust against decoherence, therefore we call their superposi- 
tions atomic Schrodinger cat states. The initial fast regime of 
the time evolution is associated with the process of decoher- 
ence, and it is directed towards the statistical mixture of the 
constituent coherent states of the original state for most of 
the initial conditions. However, certain superpositions, called 
symmetric, exhibit exceptionally slow decoherence. By intro- 
ducing a new measure, we generalize the usual decoherence 
scheme regarding the evolution of the state to account also 
for the symmetric case. To stress the fact that the environ- 
ment preserves symmetric superpositions much longer than 
the other ones, we present Wigner function images of the de- 
coherence of a suitably oriented four component cat state. 

PACS: 03.65.Bz, 03.67.L, 42. 50. Fx 



I. INTRODUCTION 

The apparent lack of a superposition of macroscopically 
distinct quantum states (Schrodinger cats) has been an in- 
teresting and vivid problem since Schrodinger's famous paper 
[jjj . A successful approach, initiated by Zeh Q and developed 
by Zurek obtains the loss of quantum coherence as the 
consequence of the unavoidable interaction with the environ- 
ment. Theoretical studies in this framework have investigated 
a variety of model systems usually coupled to a collection of 
harmonic oscillators as an environment. Fundamental work 
has been done on this subject in fi|-|ic|], for reviews see JulH . 
Important experiments have been carried out during the last 
years ]l3]-|l5| and the possible production of long-lived quan- 
tum superpositions has gained wide attention because of their 
expected use in quantum computation Jlij ]. 

In the present paper we investigate a system which is a can- 
didate for the experimental study of decoherence and possibly 
also for practical applications. The model to be described in 
detail in section [n], consists of several identical two-level atoms 
(the system) interacting with a large number of photon modes 
in a thermal state (the environment). It has the advantage 
that it is simple to make the correct transition from a micro- 
scopic system to a macroscopic one by increasing the number 
of atoms. We analyze the evolution of the reduced density 
matrix of the atomic system by the appropriate master equa- 
tion Jl7| — p~9|j - One of the main results of the present work is 
that we show how one and the same solution of this master 
equation describes both decoherence and dissipation. 

By analytical short-time calculations we show in section 



[II that the atomic coherent states j20j of our system are ro- 
bust against decoherence caused by the realistic interaction 
we consider. Since this behavior justifies that the superpo- 
sitions of atomic coherent states are relevant with respect to 
the original problem of Schrodinger, such a superposition is 
rightly called an atomic Schrodinger cat state [|2l|— 123|| . We 
also note that there are several proposals for the experimen- 
tal preparation of these type of states p4|j2^ ]. 

Next, in section ^ we present the decoherence and dissi- 
pation properties of atomic Schrodinger cat states, based on 
numerical computations of their time evolution. The time 
scales of decoherence and dissipation differ by orders of mag- 
nitude. Using this fact, we show how one can make a clear 
distinction between these two processes despite of the inter- 
play between them, and we define the characteristic time of 
decoherence. This decoherence time strongly depends on the 
initial conditions, notably, it is particularly large for a special 
set of initial cat states |^2|,^6j. This will be termed as slow 
decoherence in contrast with the general case which will be 
referred to as rapid decoherence. 

The interplay between decoherence and energy dissipation 
is the most appreciable in connection with the concept of 
pointer states M. It will be shown that when the decoher- 
ence is rapid, then the constituent coherent states of the 
initial state are pointer states to a very good approxima- 
tion. However, when there is enough time for dissipation, 
i.e. when decoherence is slow, then the initial atomic coher- 
ent states themselves evolve into mixtures, and therefore a 
refined scheme of decoherence holds. 

In order to underline the contrast between rapid and slow 
decoherence we superpose four atomic coherent states corre- 
sponding to the vertices of a suitably oriented tetrahedron. 
The time evolution of this four component cat state will be 
studied by the aid of the spherical Wigner function in section 
|y|. As it is expected, the interaction with the environment se- 
lects that pair from the initial superposition which constitutes 
a long-lived cat state. 



II. DESCRIPTION OF THE MODEL 

We consider a system of identical two-level atoms interact- 
ing with the environment of macroscopic number of photon 
modes. With dipole interaction and in the rotating wave ap- 
proximation the total system is described by the following 
model Hamiltonian: 

H = TiuJnJz + hu) k a\a k + ^9k ( a tJ- + a kJ+) , (1) 

k k 

where w a is the transition frequency between the two atomic 
energy levels, u>k denote the frequencies of the modes of the 
environment and g k are coupling constants. J+ , J_ and J z are 
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dimensionless collective atomic operators obeying the usual 
angular momentum commutation relations |E7|. This Hamil- 
tonian leads to the interaction picture master equation for the 
reduced density matrix of the atomic system, p Jl7],[l8|: 

^ = ~\ i(n) + 1) (J+ J- P (t) + p(t)J + J- - 2J^p(t)J+) 
-\ (n) (J-J+p(t) + p(t)J~J+ -2J+p(t)J-). (2) 

Here (n) is the mean number of photons in the environment 
and 7 denotes the damping rate. Note that the same mas- 
ter equation can be obtained by considering a low-Q cavity 
containing Rydberg atoms jigj ] . 

If the state of the atomic system was initially invariant with 
respect to the permutations of the atoms, i.e. it was a super- 
position of the totally symmetric Dicke states |^7| , the dipole 
interaction described by Hi nt — ^ fe figk (a\ J- + auJ+) in (|l]) 
would not destroy this symmetry. Therefore we may restrict 
our investigation of the N atom system to the totally symmet- 
ric N + 1 dimensional subspace of the whole Hilbert-space of 
the atomic system. This subspace is isomorphic to an angular 
momentum eigensubspace labeled by j = N/2. This model 
has been proven to be valid in cavity QED experiments with 
many atoms, as reviewed in p"s| ]. 

The environment as a static reservoir (represented by 
the thermal photon modes) continuously interacts with the 
atomic system influencing its dynamics. As it is obvious, the 
dissipation of the energy leads to thermal equilibrium in the 
system, corresponding to the stationary solution of the master 
equation (|^). However, as it will be shown here, the same mas- 
ter equation describes also a much more interesting process. 
The continuous "monitoring" || of the atomic system by the 
environment results in the total loss of the coherence of the 
quantum superpositions in the system. This decoherence pro- 
cess is generally extremely fast compared to the dissipation, 
except for specia l initial conditions which will be discussed in 
subsection IV A. 
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III. 



THE INITIAL STAGE OF THE TIME 
EVOLUTION 



In this section we apply general concepts to our system 
in order to find the initial states for the master equation (□) 
which are relevant to the original problem of Schrodinger M 
concerning the unobservability of macroscopic superpositions. 

First we consider the short-time behavior of the total sys- 
tem. The (by assumption pure) system+environment state, 
\9) can be written in the Schmidt representation p8|- po[|l2] ] 
at any time as 



!*(*)> = Yl v^(*) !*>*(*)> IM*)>, 



(3) 



where \ipk) and are elements of certain orthonormal bases 
(Schmidt bases) of the system and the environment, respec- 
tively. 

At zero temperature the photon field of the present model 
is in its pure vacuum state |0), therefore the initial state fac- 
torizes as 



i.e. po(0) = 1 and Pk(0) = for k ^ 0. Due to the interaction 
this product state evolves into a more general Schmidt sum 
like Eq. ^ , or in other words it turns into an entangled state. 
The time scale of this entanglement or de-separation can be 
obtained |29[ in leading order of time as 



Po(t) 



■At* 



with 



A= J2 \(M0)\(M0)\H int \M0))\0)\ 



(5) 



(6) 



This quantity can be called the rate of entanglement. Us- 
ing the explicit form of the interaction Hamiltonian H lnt , a 
straightforward calculation leads to 



A = C(J+,J_) := - (J+) (J-) . 



(7) 



i.e. in our system the rate of entanglement is the normally 
ordered correlation function of the operators J + and J_ . 

Let us turn to the case of finite temperatures, when the 
total system has to be represented by a mixed state even at 
t = 0. The linear entropy, defined as 



(8) 



Sun = Tr(p - p ), 



can be regarded as a relevant measure of decoherence ]12| , pl| . 
Restricting ourselves again to the initial regime of the time 
evolution, we can make use of the master equation ^ and 
calculate the time derivative of the linear entropy at t = 0: 

(nr) t=0 = 7 ((n> c {J ~' J+) + (<n) + 1)c (J+ ' J - }) • (9) 

The normally (antinormally) ordered correlation func- 
tion, C(J+,J~) (C(J-,J+)), disappears in the eigenstate 
\j, m — —j) m — j)) of J- ( J+). However, the collective 
atomic operators J_ and J+ have no simultaneous eigenstates 
which would annullate the right hand side of Eq. (H). Nev- 
ertheless, we are going to show that if the number of atoms 
N — j/2 is large enough, then the correlation functions in Eq. 
([]) are negligible in a class of states called atomic coherent 
states |Q. These states are labeled by a complex parameter 
r = tan(/3/2) exp(— i<j)) (for the angles j3 and (f> see Fig. [I]) and 
can be expanded in terms of the eigenstates of the operator 
J z (Dicke states) as 



2j 
j + m 



j+m 



(l + |rP)i 



(10) 



For large j, the atomic coherent states are approximate 
eigenstates of the operators J_ and J+ This state- 

ment is understood in the sense that the square of the cosine 
of the angle a between |r) and J_ \t): 



\{r\J-W)\ 



(r|r) (r\J+J-\r) 



(11) 



differs from unity by a factor which scales as (jr 2 ) -1 . Thus 
a becomes negligible in the j — » oo limit for finite r | p6[ . 
The same statement holds for the operator J+, therefore both 
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correlation functions in Eq. (^j) are indeed negligible in the 
atomic coherent states 

This suggests that the atomic coherent states are rather 
stable against the decoherence induced by the photon modes, 
i.e. they can serve as a model of classical-like macroscopic 
quantum states. This result is in analogy with the stability 
of the oscillator coherent states obtained in |3^] . 

Two such states, |ti) and \t2) can be considered as macro- 
scopically distinct, whenever the distance between the param- 
eters ti and T2 is sufficiently large on the complex plane. This 
implies that the coherent superposition of these states, which 
will be introduced in the next section, yields an appropriate 
model of the original paradox of Schrodinger. 



IV. DECOHERENCE OF ATOMIC 
SCHRODINGER CATS 

Based on the results of the previous section, the superpo- 
sitions 



1*12} = 



+ M 



a/2(1 + Re <ri|r 2 )) 



(12) 



will be called atomic Schrodinger cat states pl| , p2[|3^ ] , see Fig. 
|lj Now we are going to present our results on the decoherence 
and dissipation dynamics of these type of states. 



A. Time scales 

A typical result of the numerical integration of Eq. (Q) is 
that the time evolution of the states given by Eq. (Q) can 
be characterized by two different time scales, as illustrated by 
Fig. ^, where the linear entropy and the energy of the atomic 
system is plotted versus time. 
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FIG. 2. The two regimes of the time evolution. (Initially: 
ri = tan7r/4, T2 = 0.) The number of atoms is N = 500 and 
the average number of photons is (n) = 1, id ~ 6 x 10~ 5 /7 
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FIG. 1. Scheme of an atomic Schrodinger cat state defined 
by Eq. (p^). The points labeled by n and T2 represent the 
corresponding atomic coherent states on the surface of the 
Bloch-sphere. The angles defining the r parameters are also 
shown. 



As we can see.there exists a time instant t& (marked with 
an arrow in Fig. Eh when the character of the physical process 
changes radically. Initially Sunit) increases rapidly while the 
dissipated energy of the atoms is just a small fraction of that 
part of the energy which will eventually be transferred to 
the environment. On the other hand, for longer times t 3> 
td both curves change on the same time scale. The energy 
of the atomic system decays exponentially as a function of 
time allowing for identifying the characteristic time of the 
dissipation tdiss, with the inverse of the exponent. (We note 
that in Fig. ^| the plotted time interval is much shorter than 
tdiss, thus the exponential behavior is not seen.) More detailed 
calculations have shown that for high temperatures energy 
and linear entropy exhibit similar exponential behavior in the 
second regime of the time evolution. Their exponents coincide 
with 2-3% relative error. This implies that the initial stage 
of the time evolution is dominated by decoherence while after 
td the dissipation determines the dynamics. Accordingly we 
define the characteristic time of the decoherence as the instant 
when the slope of the curve Su n (t) decreases appreciably. We 
note that id defined in this way is in accordance with the 
decoherence time defined previously in [E^ for a specific initial 
state. 

It is remarkable that although a few hundred atoms do not 
really constitute a macroscopic system, the difference of the 
time scales is obviously seen in Fig. H. It is generally true 
that the larger N is, the more naturally and sharply the time 
evolution splits into two regimes. 

Now we turn to the investigation of the dependence of the 
decoherence time ta on the initial conditions. 
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FIG. 3. The dependence of the characteristic time of the 
decoherence on the parameters of the initial Schrodinger cat 
state: n = tan/3i/2, r 2 = tan/32/2. The number of atoms is 
N = 50 and the average number of photons is (n) = 3. 

Fig. ^ shows the contour plot of the decoherence time ver- 
sus the parameters 0i and P2 (see Fig. |l]) of the initial atomic 
Schrodinger cat state (|l2|). We have set tf>i = (p% = for sim- 
plicity. As we can see, the effect of decoherence is remarkably 
slower when /3i » @2 which was expected since in this case 
the overlap of the two initial coherent states is not negligi- 
ble, so these states can not be considered as "macroscopically 
distinct". Much more surprising is the fact that cat states 
which were initially symmetric with respect to the (x, y) plane 
(i. e. P\ ~ 7v — fa ) also decohere slower |2^|^| , but it is in ac- 
cordance with the analytical estimations of Braun et. al. Ji(]] 
In the following sections we shall refer to these states as sym- 
metric ones. 



B. The direction of the decoherence 

We saw in the previous subsection that the interplay be- 
tween decoherence and dissipation is reflected in the time evo- 
lution of the superpositions given by Eq. ([l2]). In this sub- 
section we shall focus on the direction of the process resulting 
from the dynamics governed by the master equation (Q) . 

According to the general theory of environment induced de- 
coherence ]3] Jri| , ^2|j31p3| , the interaction with a large number 
of degrees of freedom selects naturally the so-called pointer 
basis Q in the Hilbert-space of the system subject to decoher- 
ence. This process is claimed to favor the constituent states 
of the pointer basis in the sense that the system is driven 
towards a classical statistical mixture of these states. Thus, 
from the present point of view p(id) is the relevant quantity 
to be examined. 

Recalling the analytical results of sec. frill, it seems plausible 
to expect that the atomic coherent states (|10[) will be pointer 
states. 

By introducing the density matrix which corresponds to 
the classical statistical mixture of the initial coherent states: 



the expected scheme of the decoherence reads: 

|*12)<*12| ->Pci(ti,t 2 ). (14) 
We shall refer to p c \ as the classical density matrix. 
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FIG. 4. The linear entropy Sn n and the distance D between 
p and p c i (defined by Eq. (^|)) in the case of a rapidly de- 
cohering Schrodinger cat state (ti = tan7r/4, t% = 0). The 
number of atoms is N = 500 and the average number of pho- 
tons is (n) = 1. 

The distance between the actual density matrix p(t) and 
p c i, defined with 



Z?(p(t),p c i)=Tr[(p(i)-p cl ) 2 



(15) 



is always decreasing fast. Except for the case of slowly deco- 
hering cat states which will be discussed below, D(p(t),p c \) 
reaches its minimal value at the decoherence time, see Fig. 

This minimal value is very close to zero implying that the 
density matrix of the system at this instant is nearly the same 
as the classical density matrix (^) . This fact justifies the def- 
inition of the characteristic time of the decoherence in section 



IV A, and it is in excellent agreement with the decoherence 
scheme (|l4|). 

Due to the exceptionally slow decoherence, we have to mod- 
ify this picture if the initial state is a symmetric superposition. 
In this case the decoherence time is so long that the atomic 
coherent constituents of the initial state are also appreciably 
affected by the time evolution until the decoherence time, td- 
The state of the atomic system at td will be a mixture, which 
is the same as if the system had started from p c \ at t=0. In 
other words the evolution follows the modified scheme: 



|*12>(*1 



Pci(n,r 2 ,t) 



(16) 



pd(ri,r 2 ) = - (|Ti)(n| + \t2)(t 2 \) . 



(13) 



where the time dependent classical density matrix p c i(ti, T2, t) 
is the one which would evolve from the statistical mixture ( |l3| ) 
Pci(tl, T2) = p c i(Ti,T2,i = 0) according to the same master 
equation (^|) as the actual atomic density matrix. The dis- 
tance between the time dependent classical density matrix, p c i 
and p(t) becomes negligible at ta, and asymptotically reaches 
zero for long times in the case of all the initial conditions. 
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V. WIGNER FUNCTIONS OF FOUR 
COMPONENT SCHRODINGER CAT STATES 

The results of the previous section have shown that both 
the characteristic time and the direction of the decoherence 
strongly depend on the initial conditions. In this section we 
illustrate this fact by tracking the decoherence of the super- 
position of four atomic coherent states: 



i) 



In) + |r 2 ) + [73) + |t 4 ) 



(17) 



Since four points on the surface of a sphere are not distin- 
guished with respect to each other if and only if they are the 
vertices of a regular tetrahedron inscribed in the sphere, we 
set the components of 1*1234) according to this pattern. On 
the other hand, the 2 axis is distinguished in the present model 
because of the form of the Hamilton operator ([j]), therefore 
we orient the tetrahedron with one edge parallel to the z axis 
and the opposite edge parallel to the y axis, see Fig. || Al- 
though we have in principle two substantially different ways 
of considering the state represented by Fig. [] as the super- 
position of two atomic Schrodinger cat states, according to 
the results of the previous section one expects that the en- 
vironment naturally selects one of these possibilities via the 
different time evolutions: the quantum coherence between the 
components of the symmetric pair I* 12) oc |n) + |ts) disap- 
pears slowly, while all the other pairs are rapidly decohering 
superpositions. 




tion over the unit sphere (which is the appropriate phase space 
in the present case) being in a linear one-to-one correspon- 
dence with the density matrix of the atomic system. 

Following the construction of Agarwal B4], the Wigner 
function corresponding to a reduced density matrix p of the 
atomic system is defined as the expectation value of the op- 
erator kernel A(8, </)): 



with 



W{0,4>) = Tr(pA(0,4>)), 



K=0 Q=-K 



(18) 



(19) 



Here the Tkq are the multipole operators and Ykq denote 
the spherical harmonics |3^] . For previous applications of the 
spherical Wigner function see |21p^ , ^6| -^| . 

The Wigner function ( |l8| ) suggestively maps the time evo- 
lution of the state (JIt]) onto the unit sphere, as shown in Fig. 
^| Here we present the state of the atomic system at three 
time instants, and we plot the corresponding Wigner func- 
tions both as a polar plot and as a contour plot. Dark shades 
mean negative, light shades mean positive function values. 
The four positive lobes, pointing from the center to the ver- 
tices of the tetrahedron shown in Fig. ^, correspond to the 
four atomic coherent states in (|l7]). Due to the dissipation 
all these lobes will move slowly downwards. The initial inter- 
ference pattern (Figs. ^](a) and (b)) has the regularity of the 
tetrahedron, there are equally pronounced oscillations along 
all the edges, representing the quantum coherence between 
the coherent states. 

Figs, ^(c) and (d) depict the situation after a time which 
is short in the sense that the the shapes of the lobes of the 
coherent states are not appreciably affected (no dissipation), 
but the interference is already negligible between them, except 
for the single pair along the vertical edge of the tetrahedron. 
As it is seen from Fig. |E] this is the pair which represented 
initially the symmetric atomic Schrodinger cat state oc |ri) + 
T2) in (|l7|). The coherence between the components of this 
pair of states is nearly unaffected as shown by the strong 
oscillations. 

A qualitatively different stage of the time evolution is 
shown in Figs. n(e) and (f) at a later time. The coherent 
constituents are already affected by the dissipation, the up- 
permost one rather strongly, but the quantum coherence be- 
tween the components of the symmetric pair is still present. 
On the contrary, the interference between all the other com- 
ponents has already disappeared. 



FIG. 5. Phase space scheme of the 4 component cat state. 
The atomic coherent states constituting the superposition 
(0) are represented by the points labeled by n, ta. 
They are arranged to form the vertices of a tetrahedron as 
shown. 

We are going to visualize the decoherence process of 1*1234) 
by the aid of the spherical Wigner function. It is a real func- 
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FIG. 6. Wigner view of the decoherence of the 4 component cat state. We plot the spherical Wigner function ( |l8| ) both as 
a polar plot (a, c, e) and as a contour plot (b, d, f). A polar plot is obtained by measuring the absolute value of the function 
in the corresponding direction, and the resulting surface is shown in light where the values of the Wigner function are positive, 
and in dark where they are negative. Similarly, light shades of the contour plot correspond to positive function values, while 
dark shades mean negative function values. Plots (a, b) show the initial state, i.e. the superpositon IW1234) (Eq. (|l7|)) as shown 
in Fig. |B|. In (a) we also label the lobes corresponding to the initial coherent constituents. Plots (c, d) show the spherical 
Wigner function at t = 0.015/7, while plots (e, f) correspond to t = O.O4/7. 
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In view of the present results, if the initial state of the 
atomic system is a superposition of coherent states so that 
there are symmetric pairs of coherent states in the expansion, 
then the coherence between the components of these sym- 
metric pairs will survive much longer than between any other 
terms. 



VI. CONCLUSIONS 

We have investigated the decoherence of superpositions of 
macroscopically distinct quantum states in a system of two- 
level atoms embedded in the environment of thermal photon 
modes. Utilizing the Schmidt decomposition and the linear 
entropy, we have shown that atomic coherent states are robust 
against decoherence, both at zero and non-zero temperatures. 
This result is in analogy with the harmonic oscillator case 
and justifies the definition of atomic Schrodinger cat states as 
superpositions of atomic coherent states. 

By solving the master equation ([]) we have identified two 
different regimes of the time evolution with the help of the lin- 
ear entropy. The first one is dominated by decoherence while 
the second one is governed by dissipation. Based on sev- 
eral computational runs focusing on the characteristic times 
it was found that td decreases much faster than tdi as as the 
function of the number of atoms, N. Consequently td be- 
comes many orders of magnitude smaller than the charac- 
teristic time of dissipation for macroscopical samples, and 
even for e.g. N = 500 atoms and an average photon number 
(n) = 1 the ratio tdiss/id is around a few hundred depending 
on the initial conditions. However, there are very important 
exceptional cases, called slow decoherence, when the atomic 
coherent states constituting the initial atomic Schrodinger cat 
state are symmetric with respect to the equator of the Bloch 
sphere. 

Using a new measure D, we have shown that at the char- 
acteristic time of decoherence the system is always very close 
to the state described by the time dependent classical density 
matrix. Apart from the exceptional case of slow decoherence, 
the coherent states appearing in (jl2]) are approximate pointer 
states. When due to its symmetry the initial cat state is a 
long-lived superposition, also its constituent coherent states 
have time to transform into mixtures until td- We have given 
a modified scheme of decoherence which is valid also for slow 
decoherence. Accordingly, the orthonormal pointer basis con- 
sists of the eigenstates of p c i(Ti,T2,td)- 

We have demonstrated the important difference between 
rapid and slow decoherence by tracking the time evolution of 
a four component superposition with the help of the spherical 
Wigner function. The initial interference pattern having the 
symmetry of a tetrahedron rapidly disappears except for the 
single slowly decohering pair. 
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